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ABSTRACT 

Dynamic solutions combining optical, laser and Goddard Range and 
Range-Rate (GRARR) observations of GEOS-II have yielded definitive 
results for center of mass station coordinates for three laser and four 
GRARR tracking sites. Comparisons with independent solutions and in- 
ternal consistency suggest an accuracy of 10 meters or better in most 
cases. Comparisons of inter site distances on a datum suggest accuracy 
equal to or better than that usually associated with first order surveys. 
The outstanding consistency obtained with combinations of optical and 
scale-providing data suggest that the value of GM (= 3.986013 x lo^"* 
m^/sec^) obtained by JPL from deep space tracking is highly accurate. 


preceding page blank not filmed 

CONTENTS 

P age 

ABSTRACT Ui 

1.0 INTRODUCTION 1 

2.0 STATION POSITION ESTIMATION 3 

3.0 CONCLUSIONS 19 

REFERENCES 23 

TABLES 

Table Page 

1 GRARR and Laser Station Positions 2 

2 Orbital Elements of GEOS-II 4 

3 Estimated Station Positions, and Derived Positions 

Based Upon Optical Control Points. 6 

4 Arcs Used In Solutions 14 

5 Estimated Carnarvon Range Biases for Simultaneous 

GRARR and Laser Passes — Meters 16 

6 Estimated Carnarvon Laser Station Coordinates 17 

7 Carnarvon Station Positions in an Inter station Constrained 

Optical/GRARR/Laser Solution 18 

8 Differences in the Intersite Distances for First Order 
Surveyed Positions and Electronic Instruments on the 

NAD (meters) [Survey-Satellite] 20 

9 Comparison of Satellite- Derived and Surveyed Intersite 

Distances for Australia (meters) 21 


V 


FIGURES 


Figure Page 

1.1 Subsatellite Plots For ROSRAN 7 

1.2 Subsatellite Plots For MADGAR 8 

1.3 Subsatellite Plots For CARVON 9 

1.4 Subsatellite Plots For ULASKR 10 

1.5 Subsatellite Plots For GODLAS 11 

1.6 Subsatellite Plots For WALLAS 12 

1.7 Subsatellite Plots For CRMLAS 13 


Vi 


GODDARD RANGE AND RANGE RATE AND LASER 
STATION COORDINATES FROM GEOS-R DATA 

1.0 INTRODUCTION 

Definitive center of mass coordinates for the Tananarive, Mads^ascar; Rosman, 
North Carolina; Fairbanks, Alaska; and Carnarvon^Australia Range and Range- 
Rate sites, and the Wallops Island, Virginia; Goddard Space Flight Center (GSFC) 
and Carnarvon Laser sites have been estimated from GEOS-II data. Comparison 
[ 1] with our previously obtained optical coordinates indicates a high degree of 
consistency. Laser range, GRARR range and range rate, and optical data are 
all consistent; the station positions obtained agree with the optically derived 
results to 10 meters or better. The GSFC laser position also ^rees well with 
recent SAO estimates [2, 3] . For an example of the precision obtained, the 
surveyed interstation distance between the Fairbanks GRARR and the Edinburg, 
Texas MOTS 40 camera agreed with the dynamically obtained value to about 
nine meters. The chord length is almost 5.4 million meters. 

Twenty-nine two -day arcs containing a total of about 11,500 optical observations, 
and 70 laser and 180 GRARR range and range rate passes were used in the solu- 
tions. The data period ran from March, 1968 to May, 1969. Solutions were made 
usii^ differing weights for the various data types to evaluate the sensitivity of 
the solutions to data distribution and to reveal inconsistencies. 

The estimated coordinates are presented in Table 1. Of special interest are the 
new center of mass station coordinates for sites in western Australia. The 
relationship of local to center of mass coordinates can now be used to convert 
the survey coordinates of the Unified S-Band and C-Band Radar Systems at 
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Table 1 

GRAKR and Laser Station Positions 



(Coordinates referred to an ellipsoid of flat ening 1/298.255 and semi>roajor axis equal to 6,378,155 meters. 
These parameters are consistently used throughout the text.) 
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Carnarvon to the center of mass system. Prior to our analyses accurate center 
of mass coordinates were not available for these stations. 

Both the 1969 SAG American Geophysical Union (AGU) (2) and 1969 SAG Standard 
Karth (SF.) (3| Gravity models were used. The near identical results obtained 
indicate that our estimation scheme is relatively insensitive to gravity model 
error. 

These two gravity models produce orbits which differ along track by as much as 
40 m [4] at times, but this error was not translated into significantly differing 
station position estimates for reasons given below. 

2.0 STATION POSITION ESTIMATION 

The GEOS-II satellite orbital specifications are presented in Table 2. The orbital 
arcs used were chosen to provide the geometry needed for accurate station esti- 
mation. Resonance is a problem [4] for the GEOS satellites, but was dealt with 
by our use of multiple, short (2 days) arcs. 

Oui' approach to the station estimation problem was dynamical using Cowell’s 
method to solve the equations of motion. The solutions contained a large number 
of SAO Baker-Nunn observations in addition to the (primarily continental United 
States) STADAN and SPEOPT optical observations. Because of the large number 
of electronic compared to optical data points taken per pass, we selected elec- 
tronic data so that the total per pass was approximately 50. 

Based upon our experience [ 1] , we have concluded that the use of two-day arcs 
provides an accurate means for recovering station parameters. Such orbital arcs 
are long enough to have a large amount of data, but not so long that error growth 
becomes excessive. 
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Table 2 

Orbital Elements of GEOS-Il 


Epoch 

Apogee Height 
Perigee Height 
Eccentricity 
Inclination 
Anomalistic Period 


GEOb-II 
April 2«, 1968 
1569 Kilometers 
1077 Kilometers 
0.03 

105.8 Degi*ees 
112.1 Minutes 
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Most Daker-Nunn station positions were held fixed at their 1969 AGU values [2] . 
The STADAN and SPEOPT and certain Baker-Nunn optical stations were also 
held fixed at their GSFC 1970 values [5] . As an internal check, the Madagascar 
GRARR station coordinates were recovered with and without Tananarive MOTS 
40 :>pttcal data. Table 3 compares the present solutions with coordinates derived 
via survey data from our previous optical solutions. The agreement is very 
good, except for the cases of Alaska and Carnarvon. In these cases the closest 
optical stations were several thousand kilometers away. The new Alaska and 
Carnarvon positions are thus important for deriving positions of other stations 
in these areas for which geodetic-quality /quantity data is not available. Table 3 
also shows that our solutions were relatively unaffected by gravity model error. 
The good results were obtained by having (where possible) passes on all sides of 
the stations and in opposite directions. Figures 1.1 through 1.7 present the 
geometry of the passes in the present solutions. The length of the lines indicate 
the pass lengths. Although in the case of the Goddard and Wallops Lasers the 
coverage was poor, good results were obtained because the wealth of data from 
North American Datum (NAD) optical stations prevented large satellite position 
error over the United States. 

Table 4 presents the arcs used in these solutions and the number of observations 
and passes for each of the stations. These arcs were carefully selected in order 
to maximize the worldwide optical coverage while at the same time permitting a 
large number of electronic passes per recovered station. 

The validation of the GRARR data presented certain problems. Thanks to the 
extensive evaluation of the GRARR data by John Herbert of GSFC [6] it was felt 
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Table 3 

Estimated Station Positions, and Derived Positions 
Based Upon Optical Control Points 


Stations 

Solution 

Latitude 

K. Longitude 

Height 

(meters) 

MADGAR 

derived* 

1 

o 

o 

14'.' 8 

47“ 18' 

11'.'4 

1381 

(GRARR) 

GSFC dynamic 

(AGU)2 

14'.'5 


ll'.*4 

1382 



(SE)^ 

14'.'5 


11'.’3 

1382 

ROSRAN 

derived 

35“ 11' 

45 '.'4 

227“ 7' 

26'.'0 

814 

(GRARR) 

GSFC dynamic 

(AGU) 

45 '.'6 


25'.' 8 

819 



(SE) 

45'.'6 


25'.’ 8 

819 

ULASKR 

derived 

64“ 58' 

19'.'3 

212“ 29' 

10'.'9 

333 

(GRARR) 

GSFC dynamic 

(AGU) 

19'.'0 


12'.' 1 

340 



(SE) 

19'.'0 


12'.' 1 

340 

GODLAS 

derived 

CO 

CD 

o 

o 

13'.'9 

283“ 10' 

18’.'4 

-2 

(LASER) 

GSFC dynamic 

(AGU) 

14'.' 1 


18'.'4 

3 



(SE) 

14'.' 1 


18'.'5 

4 


SAO determinations 

(AGU) 

14'.'2 


18'.'2 

5 



(SE) 

14 '.'2 


18'.'3 

9 

WALLAS 

derived 

37“ 51' 

35'.'7 

284“ 29' 

23*.' 9 

-53 

(LASER) 

GSFC dynamic 

(AGU) 

36’.'0 


24'.' 0 

-60 



(SE) 

36'.'0 


24'.' 0 

-59 

CRMLAS 

derived 

-24“ 54' 

16'.'7 

113“ 41* 

58'.*2 

2 

(LASER) 

GSFC dynamic 

(AGU) 

16'.'4 


57'.' 9 

-5 



(SE) 

16'.'4 


57'.' 8 

-3 


1. Derived using survey data for the GRARR Station and applying shift from local 
datum to center of mass system as obtained in previous optical solutions. 

2. Dynamically estimated using the SAO AGU 1969 gravity model. 


3. Dynamically estimated using the SAO 1969 Standard Earth gravity model. 
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Figure 1.4. Subsatellite Plots For ULASKR 
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Table 4 

Arcs Used In Solutions 


1B6B ARCS 



Uf^iSKH 



MAlXiAH 


MOSMAN 

WALLAS (iODLAS 

OpUfal' 

Date 

No. of Obs, 


No, of Obs. 


No. of Obs. 

No. of 

No. of 




.. No 

. of 


No 

. of 


No. of 

( }bs 

No. of Ol)S. 

No. of 

No. of 


rUnKC Passes 

llanKi* 

Pa.*... H.nK.. ' 

llange p jjgpg 

Mange 

It, mg,. 

Passes 

obs 

1 2-:! Oh 

213 MO 

3 




IHO 

1 H(» .3 

49 i 

2 


311 

4/20-27/Oh 

I3f. Ml 

A 


1 


313 

313 1 

119 

1 


.320 

5/7-h'Ok 

172 172 

2 




102 

102 1 

IhO 

2 


172 


132 132 

2 




107 

100 1 

192 

2 


351 

0^0* 10/ OH 

202 202 

3 




277 

277 4 




393 

0/Jl »U/OH 

1 H2 1 H2 

3 




150 

130 5 

;,HH 

3 



0/H-l/iAiH 

2.i2 232 

1 




190 

1 07 3 




Oil 

0/ 10-17/OH 

130 130 

2 




103 

191 3 




HI 1 

0/21-22/08 

271 271 

5 




203 

279 1 

122 

3 


75 1 

6/23-2 l/CH 

2h5 2h5 

a 




229 

229 1 




020 

0/21-25/00 



in 

112 

2 




09 

1 

132 

0/27-2H/0H 

103 103 

3 


57 

2 




3h 

1 

3 HO 

10/4-r)/OH 

200 200 

3 

105 

107 

2 




73 

1 

2h0 

U)/0-7/6H 

271 271 

1 

1h(» 

IHO 

1 






309 

10 H-9/6H 

07 07 

1 

177 

iHM 

1 




31 H 


507 

10 '2 1-22/OH 

132 132 

<) 







312 

2 

115 

JO/23-24/OH 

202 202 

3 







137 

2 

3Ui 

niTADH 

310H 2090 

IH 

570 

(133 

M 

2370 

2393 31 

259(i 

13 1217 

9 

7h21 

t 

















1960 

ARCS 







1 LASKIM 


(’AHVON 


CM MIAS 

001)1. AS* 


OpUeal* 


No. of Obs. 



No, of Obs. 


No of 


No. of 



Date 


No, 



No. ol 

Obs. 

No. of 

Obs. No. 

of 

No. of 


Hair’ 

Passt 

Mange 

Mange 

Mate 

llanuf 

Passes 

,, Passes 

Mange 

Obs. 

3/2-3/00 

127 127 


5 

15H 

I3H 

1 

99 

2 

01 

1 

150 

3/3-0/09 

121 121 


3 

02 

129 

2 

227 

3 

300 

4 

300 

3/11-12/69 

00 00 


3 

232 

232 

1 

1 00 

3 

105 

•1 

471 

3/13-14/00 

150 150 


.3 

314 

379 

0 

1 99 

3 

32 

1 

295 

3/17-Ih/OO 

100 100 


5 

232 

237 

5 

2M 

1 

231 

3 

224 

3/29-30/60 

02 02 


5 

170 

170 

5 





225 

3/31-4/1/69 

104 103 


4 

192 

200 

1 

loj 

3 



234 

l/H-0/69 

75 75 


2 

107 

170 

1 

HO 

2 

94 

3 

3 HO 

4/10-n/OO 

172 172 


5 

00 

M3 

2 

321 

1 

73 

1 

24H 

4/M 13/00 

123 123 


3 

09 

1 2a 

3 

2V1 

1 



230 

4/24-25/00 

1 00 1 09 


5 

150 

179 

3 

251 

2 

159 

1 

1.52 

3/3-0/09 

J(i3 103 


4 

210 

2 HO 

1 

210 

1 



514 

TOTALS 

10H4 10H3 


17 

2115 

2114 

43 

2235 

33 

1127 1 

10 

3H22 

^Station coordinates held fixed. 
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9862 
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186 









Laser 
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7205 
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71 
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that the Kosman, Alaska and Madagascar data sets during 1968 (solution 1) were 
relatively free of significant biases compared to the noise level, except for a 
four meter bias due to survey error to the ranging target at Rosman. The 
Carnarvon GRARR which was also tracking GEOS-II at this time was found to 
have significant hardware problems which prohibited its inclusion in our first 
solution. The later 1969 Carnarvon GRARR data set (solution 2) also seemed to 
have significant range biases when compared to simultaneous bias-free Carnarvon 
laser passes, indicating that cei’tain problems persisted within this GRARR sys- 
tem throughout the period of interest in 1969. Table 5 presents the estimated 
range biases for the Carnarvon radar. Herbert’s biases were determined by 
analysis of GRARR residuals based upon laser short arc reference orbits. Our 
long arc solutions were run with the laser and GRARR sites both held constrained 
and allowed to adjust independently. In the constrained solution both stations 
must adjust together by the same amount. 

Along with the station coordinates, range biases were also estimated simultane- 
ously with both loose (100 meters) and tight (4 meters) a priori sigmas on the 
range biases in our solutions. As can be seer, from the estimated range bias 
values, the long arc solutions gave good agreement with the bias values ‘obtained 
from short arc solutions. In the constrained solutions the agreement was five 

meters or better in thirteen out of seventeen passes of data. This is good since 

<• 

we were solving not only for a range bias at Carnarvon but also the station co- 
ordinates as well. Thus, given the variation in the biases between the independent 
and constrained solutions and the erratic nature of the magnitude of the biases, 
it was felt that the Carnarvon GRARR station coordinates could not be satisfac- 
torily recovered independently of the laser-GRARR interstation distance. 
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Table 5 

Estimated Carnarvon Range Biases for 
Simultaneous GRARR and Laser Passes — Meters 

1969 


Data 

Herbert 

Solution 

a priori 

a ~ 4m 

a priori .r 

= 100m 

Constrained 

Independent 

Constrained Independent 

April 8-9 

-10.3 

2.5 

-1.9 

2.2 

-6.6 


-12.9 

-7.3 

4.8 

-9.0 

3.0 

April 10-11 

-8.1 

-6.3 

4.9 

-7.0 

2.6 

April 24-25 

0.6 

2. .5 

6.0 

2.5 

2.5 

March 2-3 

-9.9 

-16.4 

-9.0 

-18.0 

-13.2 


-1.3 

-0.5 

11.0 

-1.3 

9.6 

March 5-6 

-3.7 

-6.0 

-6.6 

-6.4 

-10.9 

March 11-12 

0.0 

2.5 

4.9 

2.2 

0.8 


-11.4 

-6.8 

-0.2 

-8.2 

-4.1 


-15.8 

-11.2 

-8.0 

-12.5 

-12.8 

March 13-14 

-4.8 

-4.2 

-5.2 

-5.1 

-10.1 


-7.1 

-7.6 

2.8 

-8.6 

0.3 


-12.9 

-15.5 

-20.0 

-17.4 

-26.0 

March 17-18 

-1.5 

-6.0 

-1.7 

-7.0 

-6.1 


-5.0 

3.7 

11.7 

3.6 

9.3 


-4.7 

-2.8 

-1.6 

-2.9 

-5.6 


0.6 

2.6 

7.1 

2.6 

3.8 
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Table 6 

Estimated Carnarvon Laser Station Coordinates 


Solution Number 

Latitude <N) 

E. Longitude 

Ellipsoidal Ilei^t 
(m) 

1 

-24“ 54’ 16’.’4 

113“ 42’ 57’.’8 

-5. 

2 

16’.'5 

57'.’9 

-4. 

3 

16V4 

57'.’9 

-5. 

4 

16’.'3 

57'.'8 

-5. 

5 

16’.'5 

58'.'0 

-7. 
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Therefore the final estimated laser-GRARR station coordinates were obtained 
constraining the laser and GRARR interstation distance. The final Carnarvon 
laser /GRARR coordinates selected are very satisfactory as can be seen in 
Table 6. This table gives the laser position: 

1 - when estimated independently of the GRARR position, 

2 — when estimated from an independent optical/laser-only solution, 

3 - when estimated constrained to the GRARR position and solving for a 

GRARR range bias, 

4 ~ when estimated constrained to the GRARR site using rai^e rate only, 

and finally 

5 — estimated from laser /optical solution with the laser weight being 1/25 

of the nominal value of (6 meter)’ 
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The weight normally given to the laser range data was (6 meters)' ^ . Of course 
the laser data is much more accurate than this. We obtain fits of 2-3 meters 
(rms) on a routine basis for orbital arcs of 1 to 3 days. The large nominal laser 
weight was chosen to ensure si^’nlfli^ant weight for the optical data. As seen in 
Table 6, the agreement of solutions is highly satisfactory. Due to the GRAPH 
range biases at Carnarvon previously noted, solution 3 was adopted. The sur- 
veyed distance between the Carnarvon laser and GRARR sites should be of ex- 
tremely high accuracy since they are only about 150 meters apart. Thus, we did 
not feel that this constraint adversely affected the accuracy of our recovered 
positions. Table 7 gives the final positions adopted. 


Table 7 

Carnarvon Station Positions in an Inter station Constrained 
Optical/GRARR/Laser Solution 


Station 

Latitude (N) 

E. Longitude 

Ellipsoidal Height 
(ro) 

Carnarvon 

(LASER) 

-24“ 54' 16'.'4 

113“ 42' 57'.'9 

-5. 

Carnarvon 

(GRARR) 

-24“ 54' 11V4 

113“ 42' 58'.'9 

1. 

Carnarvon 

(Optical) 

-24. 58' 23'.’4 

113“ 43' 15'.'6 

-14. 
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Survey data have also been used to provide an important check on accuracy. > 
Tables 8 and 9 present the difference in the surveyed and dynamically recc>vered 
interstation chord distances for our dynamically recovered positions on the North 
American and Australian National Datums. The accuracy of the local surveys is 
also given. Most of the optical sites on the NAD are of first order accuracy. 

On the North American Datum the GRARR and laser intersite distances are in 
good agreement with the results for previously derived optical stations. The 
chord length differences in Table 8 indicate that the satellite derived chords are 
longer than the local survey chords by approximately a few parts per million. The 
chords to Alaska indicate approximately the same scale change but in the opposite 
direction. This is not unexpected since the Alaska Station is on a third order 
survey and ties to the or igin of the North American Datum (Meades Ranch, 

Kansas) are possibly in error. It is interesting to note that in Table 9 for the 
Australia*! Datum, the satellite derived chords are also greater in length than the 
survey chords. 

3.0 CONCLUSIONS 

It is apparent that our new GRARR and Laser station positions are of comparable 
accuracy to our previous optical station coordinate solutions. The close agree- 
ment of surveyed and satellite-derived intersite distances on the NAD indicates 
the strength and insensitivity to model errors of our method of station recovery. 

A consistent agreement with surveyed values of well under ten meters, and ten 
to fifteen meter agreement to Fairbanks from the Continental United States (with 
chords almost equalling an earth radius), has been obtained. Thus the GRARR 


Table 8 

Differences in the Inter-site Distances for First Order 
Surveyed Positions and Electronic Instruments on the NAD (meters) 
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Table 9 

Comparison of Satellite- Derived and Surveyed 
Intersite Distances for Australia (meters) 
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system is a geodetic-quality instrument. The overall RMS of fit in our solutions 
was 10 m for GRARR range and 7 cm/sec for range rate. 

It is felt that this work can only be substantially improved upon if there is further 
improvement in gravity models. But perhaps as important, we have shown that 
satellite geodesy can provide accuracy comparable to that of classical methods. 
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